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Abstract 

At submillimeter wavelengths typical gap discontinuities in microstrip, CPW lines or at 
antenna terminals, which might contain diodes or active elements, cannot be viewed as 
simple quasi statically evaluated lumped elements. Planar Schotky diodes at 2. 5 THz for 
example have a footprint that is comparable to a wavelength. Thus, apart from modelling 
the diodes themselves, the connection with their exciting elements (antennas or 
microstrip) gives rise to parasitics. Full wave or strictly numeric approaches can be used 
to account for these parasitics but at the expense of generality of the solution and the 
CPU time of the calculation. In this paper an equivalent network is derived that 
accurately accounts for large gap discontinuities (with respect to a wavelength) without 
suffering from the limitations of available numeric techniques. 

I. INTRODUCTION 

Planar circuits utilizing microstrip or coplanar waveguide are often coupled to active 
devices or antennas through electrically small gaps. However, at THz frequencies, 
practical fabrication scales limit devices, such as planar Schottky diodes [1], to dimensions 
which are comparable to a wavelength. As a consequence, the effect of the gap dimensions 
on the coupling coefficient or resonant frequency of the associated transmission line or 
antenna cannot be ignored. Nor can it be modelled as a quasi statically derived [2 and 
references there cited] single lumped element. In this paper the authors present a new 
formulation for analyzing the effect of finite gap excitations on microstrip lines. The 
technique involves separating the excited modes into reactive elements which are confined 
near to the gap and propagating quasi TEM components, which are coupled to the 
surrounding transmission line and load (antenna). The formulation is very general and can 
be applied to a large class of problems involving discontinuities in transmisssion line 
circuits. Verification of the new method has been obtained using traditional but much 
more time consuming method of moments (MoM) techniques [3, 4 for example]. 

Recent investigations [5, 6, 7] present a direct integration method to derive the spectral 
Green’s function for microstrip and strip-line. Most notably they furnish a detailed analysis 
of the singularities of this Green’s function, while enphasizing on the excitation 
mechanisms of the propagating and leaky modes compatible with strip line and microstrip. 
This technique can be used to great advantage in the modeling of sub-mm wave circuits. 
The same formalism was adopted by one of the authors, [8, 9] to derive, analytically in 
that case, the Green’s function of a slot printed at the interface between two different 
homogeneous dielectrics, and in [10] to derive ad hoc basis functions in the MoM analysis 
of coplanar waveguide excited slot antennas. 

As another application of this methodology, we investigate in this paper the input 
impedance of gap excited microstrips coupled to planar antennas. A rigorous equivalent 
network is derived that accurately accounts for the field in-and-near large-or-small (with 
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respect to a wavelength) gap discontinuities, without suffering from the lack of physical 
insight of purely numerically derived solutions. On the contrary the network approach will 
provide the guidelines on the one hand for the efficient coupling to and from the gap, and 
on the other hand for the separate characterization of the reactive energy associated with 
the feed. Full wave tools, such as Finite Difference in Time Domain or Finite Element 
Method may be necessary to characterize the minute details of a complicated diode 
structure but are typically not needed for the remaining part of the circuit. 

In our technique the spectral kernel of the Greens function is first used for the 
determination of the bound mode wave number. Then the evaluation of the residue 
contributions associated with this pole allows the derivation of an equivalent transformer 
that accounts for the gap launching coefficient on the equivalent transmission line of the 
bound TEM mode in the microstrip. By subtracting the current associated with this main 
mode on the infinite line from the total current, a source attached current is defined. This 
current, for the circuit substrates used in many sub millimeter planar devices, is mostly 
confined in a region very close to the gap. From this current a parallel gap impedance is 
derived that accurately represents the reactive effects of the gap. The same formalism 
adopted here for the microstrip configuration, can be used for characterizing a wide 
variety of other gap excited planar transmission lines and active antennas. 


II. FORMULATION: ELECTRIC CURRENTS ON THE MICROSTRIP 

The problem under investigation is shown in Fig. 1 . 
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Fig.l Infinitely extended microstrip excited with by an impressed electric field distributed over a gap of 
dimensions (5, wj. 


An infinitely extended (along x) microstrip, printed on a grounded dielectric slab (e r j, 
t r 2 =l) is considered. The microstrip cross section (w, along y) is uniform and much 
smaller than the wavelength. The structure is assumed to be excited by an impressed 
electric field e f (x,y)=ef(x,y)i x . The electric current distribution is assumed to be 
characterized by a separable space-dependence with respect to x and y; i.e., 

j(x,y)=i(x) ± . I— (1) 

V 1 v w / 

where the transverse y-dependence is chosen to verify the quasi-static edge-singularities 
and the the normalization constant, has been chosen in such a way that i(x) represents 
a net current flowing in the microstrip at any cross section of the slot line. The method 
adopted to derive the current of a gap excited microstrip has been presented in [5]-[7]. A 
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brief review with the minor variations introduced in this paper is presented in the 
appendix. As a result the electric current in the microstrip is expressed as: 


i( x )= h f 


-E> (k x )e’ jkxX 
D(k x ) 


dk x 


( 2 ) 


where E f ( k x ) is the Fourier Transform (FT) of the impressed electric field, with respect to 
the variable x, and 


OO 

D ( k x) = ik I G xx(kx.ky) J o(\wsky )dky 

-OO 


(V 


with J 0 is the Bessel function of zero order and ky) the FT, with respect to the 

variables x and y of the electric field radiated by a short electric dipole lying at the 
interface between the slab and the homogeneous half space for z<0, in the absence of the 
metallic strip. 

An interesting aspect of (2) is that the impact of the feed is all contained in the numerator 
and in most practical cases it is a simple anayitic function. When the distribution is uniform 
over the feed region (x)=-l/6 for x € (-6/2, 6/2) the excitation is said to be a 5-gap and 
E f (k x )= - sinc(k x 6/2). This case will be considered specifically in the following. The 
denominator, represents the micro strip spectral Green’s function and for a detailed 
discussion on the nature of the branch singularities of D(k x ) the reader is referred to [5, 6]. 



Fig, 2 Complex k x plane topology for thin dielectric substrates. Only the TMO pole of the grounded slab 
is assumed to be above cut off. 


The outer integral in (2) can be performed by adopting a deformation in the complex plane 
below the real axis for Re (k x )<0 and above the real axis for Re (k : x )>0 ( original path in 
Fig. 2), in order to avoid the branch cuts arising from the branch points implicit in D(k x ). 
Here it is sufficient to know that these branches are square root type or logarithmic when 
they originate from poles or from the branch (ko) of the layered Green’s function 
G£(k x ,k y ) respectively. For a thin dielectric layer configuration, (that is only the TMO pole 


529 


12 th International Symposium on Space Terahertz Technology > 


is above cut off in the grounded slab) the complex plane topology of the microstrip’s 
Green’s function, is well represented in Fig. 2 which shows the branch cuts arising 
from kj-MOi kj£j and ko. 

11.1 Case study 

It is instructive to perform a parametric study of the currents as a function of the gap 
dimensions. The information derived will lead to the current representation and equivalent 
network of the following sections. Fig. 3. shows the real part (3a) and imaginary part (3b) 
of the current due to a 1 volt uniform excitation for different gap dimensions (6) from 2 to 
18 n m. The width of the microstrip (w) and the height of the slab (h) are 8 pm and 3 pm 
respectively (h—Xg/40 at 2.5 THz). The relative dielectric constant considered is 
€ rl = 12.85. These parameters and the maximum gap dimension are typical of thin GaAs 
membrane diode technology [1] for which this analysis is derived. 




a) b) 

Fig. 3 Current on an example microstrip structure at submillimeter wave frequencies (2.5 THz) obtained 
using the present direct integration method (equation 2). Real (a) and Imaginary (b) parts are presented 
for different dimensions of the excitation gap (6 varying from 2 to 18 pm). 


Two aspects of the problem should be noted 

0 A launching efficiency issue stands out, since the amplitude of the oscillation of the real 
and imaginary part of the current at a certain distance from the source is regulated by the 
gap dimensions. As 6 tends to zero, the amplitude of the current is a maximum and tends 
to one half of the characteristic admittance of the microstrip line. As 6 grows the value of 
the current diminishes. 

ii) The reactive energy in the feed zone is strongly affected by the width of the source. 
When the gap dimensions shrink to zero, the imaginary part of the current tends to grow 
logarithmically (Fig. 3b) indicating a capacitive reactance. When the dimensions of the gap 
are large, the imaginary part of the current can be less then zero, indicating an inductive 
reactance. 

II.2 Modal and source attached currents 

The dispersion equation, D(k x )=0, defines the position of the pole of the guiding structure. 
The localization of the zeroes of D(k x ) as a function of the parameters of the structure (w, 
e r] ), can be performed numerically by means of a descent along the gradient algorithm. As 
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a starting point one can use the analytical formulas [1], knowing that the actual 
propagation constant is not drastically different. Only one pole, indicated here after as 
±k xp , is found for the thin dielectric layers we are interested in at 2.5 THz). In 

the absence of losses, the pole of the structure is found on the real axis of the k x plane (see 
Fig. 2). Knowingn its location, the integral in (2) can be approximated asymptotically for 
| x | >8/2, by the residues of the poles in k x = ±k xp . These residue contributions are 
expressed as 

i res (x)= - 2nj(sinc (k x0 - 2 )) R e~ Jk ^R=Res (^h) k - kxp (4) 

It is apparent that this current is the one associated with the quasi TEM mode that 
propagates in the microstrip. A comparison between the residue evaluation (4) of the 
integral in (2) and the original numerical integration is presented in Fig.4a and Fig. 4b, for 
two unit 5-gap excitation configurations. Outside the source region (gray zone in Fig. 4) 
the full integration of (2) and the sole residue contribution are almost superimposed due to 
the drawing scale. A source attached current can then be defined and will be indicated by 

i s (x)=i(x)-i res (x) (5) 

The source attached currents represent an overall current that contains all the separate 
branch conntributions, without calculating them separately as in [7], In this respect it is 
useful to observe that i s can be calculated for | x \ >8/2, deforming rhe original integration 
path on a path Pj (Fig. 2 where only the spectrum for Re(fcj>0 is shown ) that surrounds 
all the branch cuts on the top Rieman sheet. This is possible because in thin dielectric 
cases, the only significant singularity encountered in the deformation is the pole in ± k xp . 




Fig, 4 Comparison between asymptotic evaluation (residue only) and direct numerical integration of 
the current for a sample microstrip structure at submillimeter wave frequencies (2.5 THz). Also plotted is 
the source attached field (real and imaginary parts in 4.b and only the imaginary part in 4. a since the 
real part is almost zero and cannot be distinguished in the drawing scale). The width of the microstrip 
and the dimensions of the excitation gap (S) are A 0 /20 (6pm), \ 0 /60 (2pm) respectively. The height of 
the slab is \ 0 /40 (3pm) in Fig. 4a and \ 0 /20 (6pm) in Fig. 4b. The relative dielectric constant considered 

is (€ r j =12.85). 
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The integrand decays exponentially in the deformed path P d . Thus the numerical 
calculation of i s is much more efficient than along the original path. The definition of the 
source attached current is also analytically extended to the source region, | x | <6/2 (but 
not this way of calculating it, which must be done in accordance with the definition in 
(5)). The source attached field is also plotted in Fig. 4. In [7] it is demonstrated that each 
of the branch contributions to the current decays as x~ 3a , thus also f(x) decays this way 
or faster. It’s value in the source zone depends on the thickness of the dielectric 
considered. In the very thin dielectric slab ( h=\ o /40 ) considered in Fig. 4a the real part of 
P (x) is not noticeable on the drawing scale and thus is not reported. Also the imaginary 
part of i s (x) is very concentrated in the surrounding of the source. In Fig. 4b the same 
curves are reported for a thicker dielectric case (h=\ 0 /20). In this case the real part oii s (x) 
is noticeable, and significant for larger distances from the source. It is worth repeating at 
this point that (h=\ 0 /40) corresponds to a 3^m thick dielectric at 2.5 THz, the design 
target for this work. 

m. INPUT ADMITTANCE AND EQUIVALENT NETWORK 

The current representation in terms of source attached and mode contributions is 
particularly significant when only one bound mode is propagating in the microstrip. In this 
case it is possible to derive a rigorous equivalent network, based on a single transmission 
line, that characterizes the 6-gap excited infinite microstrip input impedance. We will see 
in section IV, that finite (loaded) microstrips can also be studied by resorting to the same 
simple equivalent network. The input admittance of a 6-gap excited micro strip can be 
expressed in the spectral domain as: 

oo 

>V-/ ^ f 1 dk, (6) 

-oo 

Accordingly, the total complex power radiated by the source can be expressed as: 

Ptot~2 I Vg | 2 yin/ (7) 

The knowledge of the microstrip characteristic impedance zo [3], together with the 
amplitude of the traveling wave in the microstrip from (4), consents to write the power 
launched in the microstrip quasi TEM mode as: 

Pm=\V g \ 2 (R + ) 2 z 0 (8) 

where R + =- 2nj( sine (k x0 ~ 2 ))R. This power actually leaves the excitation zone in form 
of a couple of outgoing travelling waves ( ± k xp ). As we will see later on, this power is the 
only one that might give rise to interactions with components of the circuit which are 
located at a significant distance, in terms of a wavelength, from the feeding zone. In order 
to derive an equivalent network that rigorously takes into account the gap effects one may 
associate all the power that is not launched into the the quasi TEM mode as being 
assigned to a lumped load indicated as z gap . The power conservation then imposes: 

Pgap-Ptol Pm (9) 
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This gap impedance is imagined as in parallel to the source and to the modal loading as in 
fig. 5a. Accordingly both P m and P gap can also be expressed as: 

P m = 2 I Pg I 2 y'm ’ Pgap ~2 I I ^.Fg<V> (10) 

thus givin the exphcit relation between the various admittance components of the network 
in fig.5a 

ygap~~yinf Tm (11) 


The expression for y m according to this definition equating the two expression for P m in 
(8) and (10) 

y m =2(R + ) 2 z 0 (12) 



Fig. 5 Equivalent network for 6 gap excited microstrip structure, a) admittance decomposition 
(y in f=y gap +y m ); b) Interpretation of the modal impedance interms of transformer connected equivalent 
transmission line . 


So that, given (12), (11) becomes an operative definition. It is useful to perform a further 
algebraic manipulation expressing the modal admittance as 


) } m 


ft 


2z, 


(13) 


where 


n =2z 0 R + = 2z 0 (2 / KjR ) sine (kx 0 1 ) 


(14) 


The introduction of the coefficient n consents to represent the modal part of the 
impedance associated to the series of two transmission lines representing the quasi TEM 
mode. These transmission line are series fed, via a transformer of turn ratio n, by the 
parallel combination of the source voltage and the gap impedance. The fall circuit is 
shown in Fig. 5b. 

The use of the transformer is a convenient wayof separating the different field behavior 
associated with the feeding and the propagation of the transmission lines. The electric 
currents in the source are oriented parallel those in the microstrip, both along x. The 
electric fields in the source and in the transmission line are orthogonal, (along x and along 
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z respectively). The value of the transformer n for the same microstrip of Fig.3 (operated 
again at 2.5 THz) is plotted in Fig. 6. as a function of the gap’s width, 6. 



Fig. 6 Transformer for 6 gap excited microstrip structure of Fig. 3 



Fig. 7 Total and gap admittance (gap and total susceptances are equal) 

The gap and total admittance for to the same case are presented in fig. 7. It should be 
noted that the gap conductance is virtually zero for the case of the thin dielectric slab, and 
the gap susceptance, coincident by definition with the total one, goes from being 
capacitive to inductive for the investigated range of 6. 

TV. NUMERICAL EXAMPLES, APPLICABILITY AND CONSIDERATIONS 
The purpose of deriving the equivalent network is to characterize the input impedance of 
finite, as opposed to infinite, microstrips. When the bound mode currents are much more 
intense than the source attached currents, at a certain distance from the source one may 
assume that the source attached currents intrinsically satisfy Maxwell’s equations for the 
finite structure. The loading of the end-points of the microstrip will generate a reflection 
that can be characterized by a load impedance. Thus the interaction between the two outer 
loads and the source can be accounted for via a complete equivalent network with 
equivalent lumped elements z// z/^ (Fig. 8) located at the endpoints of the microstrip. 


534 




12 th International Symposium on Space Terahertz Technology 


IV. 1 Numerical examples 

Two finite length microstrips printed on a GaAs membrane (e r i=12.85) have been 
investigated using the equivalent network in fig. 8. 


/ / 



Fig. 8 Equivalent network for a 6-gap excited and loaded microstrip. 

The dielectric thicness and the width of the conductor are 3 pm and 2.5 pm respectively. 
The pertinent open end impedance has been derived as in [3]. Fig.9a shows the imaginary 
part of the input impedance for a 60pm long microstrip as a function of the frequency for 
three different gap dimensions (0.5, 1.15, and 15 pm). It is found that in the three cases 
the resonant frequency is the same and significant differences are observed at larger values 
of the impedance. Also the pertinent curves obtained via a full wave Method of Moments 
are shown for the 1.15 pm and the 1 5 pm cases. The MoM tool uses piece wise sinusoidal 
sub-domain basis functions and rectangular weighting functions in the longitidinal 
dimension and assumes a rectangular distribution in the transverse dimension. The 
agreement between the network solution and the MoM solution is excellent apart for a 
minimal (0.5 %) resonant frequency shift, probably due to the different electric current 
transverse distribution assumed in the two models. 



a) b) 

Fig. 9 Impedance of a delta gap excited microstrip as function of the frequency, a) (1=60 pm); b) 
l=40pm). Also a parametric study as a function of 6 and comparison with MoM are shown. The width of 
the microstrip and the thickness of the dielectric are 2.5 pm and 3 pm respectively. The grounded 
dielectric is made of Ga-As. 

The case in Fig. 9a (1=60 pm) is geometry corresponds to a case in which the load 
impedance (2//, i=l,2 ), when reflected the excitation zone, is nearly zero. In view of this 
the gap impedance effect is relatively minimal since it is in parallel to a very low 
impedance load. The case in Fig. 9b (l=40pm) corresponds to a very high load impedance 
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at the excitation zone, and as seen the impact of the gap becomes dramatic. The resonant 
frequency as a function of the gap width shifts very significantly for 6 varying from 2.3 /urn 
to 6.9 pm and the curves can easily be distinguished. For even smaller gaps the reso nan t 
frequency shifts out of the figure since the source attached currents tend to be s ingular as 
6 tends to zero. In Fig.9a is also plotted a MoM result, for 6=3.84 pm, which again shows 
excellent agreement with out model. Finally in Fig. 9b the result that would be achieved by 
using the equivalent network without the transformer and the gap impedance introduced in 
this paper is plotted. Not only is the resonance frequency poorly modelled but also the 
absolute values of the impedance vary widely form those of the more accurate equivalent 
network analysis. 

IV. 2 Applicability 

We will assume than when | i s (x=l/2) | the accuracy, A(x=l/2) >^, in the 

evaluation of input admittance is satisfactory. In fig. 10 | i res | / 1 j is plotted as a 
function of the frequency at an observation point located at X 0 /10 from the center of the 
feed. The microstrip’s nominal characteristic impedance is 50 ohms, the dielectric 
thickness and the width of the conductor strip are 3 pm and 2.5 pm respectively. The 
dimension of the excited gap is X 0 /40 and the dielectric is composed of silicon {t r i - 1 1.7). 
It is found that | i res | / 1 i s | is always larger than 100 for frequencies up to 3 THz. Thus, 
since i s (x) decays algebraically, x~ 3/2 , any load at a distance larger than \ 0 /l 0 can be 
characterized accurately with the present equivalent network. The same is not true for 
thicker dielectrics. Also the investigated microstrip is only weakly dispersive, as the 
normalized propagation constant k xp /k 0 calculated as explained in sect. Ill varies from 2.8 
to 2.9 m the frequency range considered in Fig. 10. 

IV. 3 Considerations 

The loads. Fig. (8), typically represent active impedances. For example if they represented 
slot antennas their impedance should be calculated when both slots are radiating in phase. 
In this case the load connected at the end points of the microstrip will present a larger real 
part. Also in this case frequency shifts as significant as in Fig. 9b, can be expected. This 
has already been experimentally verified in a similar slot circuit [1 1], If the actual source 
under investigation is not well represented by a 6-gap generator, a correction in the 
equivalent network parameters can be made by simply substituting the sine function in (6), 
(14) with the relevant FT E? (k x ), provided the more realistic function e? (x) is still an even 
function, with respect to x. If ef (x) also presents some odd components, then a s ing le 
impedance should not be used, and an even and odd impedance characterization should be 
adopted. The extra degree of freedom (odd mode) can still benefit from an equivalent 
network, that can be derived with the same methodology used here for the even mode. At 
this point e f (x) could for example the result of the investigation of complicated source in 
isolation. Indeed the procedures presented in this paper are found to be particularly useful 
in connecting an isolated load or source to the remaining part of the circuit. Full wave 
tools such as Finite Difference Time Domain or Finite Element Method may be necessary 
to caharacterize the tiny details of a complicated source (like a diode) but are often not 
necessary for analyzing the remaining part of the circuit. 
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Fig. 10 Ratio between the amplitudes of the residue contribution to the currents and the source attached 
currents. The graph can be read as 1/(A(1=0.1X 0 )). The microstrip impedance is nominally 50 ohm. The 
width of the microstrip and the thickness of the dielectric are 2.5 pm and 3 pm respectively. The 
grounded dielectric is made of silicon (€ r j —11. 7). The width of the gap is Xq/40. 


V. CONCLUSIONS 

An equivalent network has been presented for the wide band characterization of the gap in 
a series fed and loaded microstrip line. The network consists of a transformer and a gap 
impedance. The transformer accounts for the launching coefficient of the mode on the 
equivalent transmission line representing the microstrip and is evaluated analytically by 
simple interpretation of the spectral expression of the microstrip Green’s function. The gap 
impedance accounts for the current distribution in the region surrounding the feed, and is 
calculated by simple spectral domain numeric integration (less then 3 minutes on a 300 
MHz Pentium II). The procedure, validated via MoM, will work for any kind of dielectric 
stratification, but is meaningful when one bound mode is above cut off in the microstrip. 
Submillimeter’s wave problems in which this equivalent network can be applied include the 
design of receivers (i.e. diodes as mixers or frequency multipliers) and generators (i.e. 
photomixer based local oscillators). 


APPENDIX: SOLUTION OF THE INTEGRAL EQUATION 

According to the equivalence principle an Electric Field Integral Equation (EFIE) is set up 
and enforced, via a point matching like procedure only on the longitudinal axis of the 
micro strip (y=0). 


If gxx (x-x 0-y')j x (x y’Jdx dy ’= - e ! (x, 0) (A0) 

E 

where £ is the entire microstrip domain. The point matching (or razor blade) as testing 
procedure is selected in order to mantain the same formalism used [8-10] for the 
characterization of slot antennas printed between two different homogeneous dielectrics. 
All the quantities are scalar, because of the x-polarisation of all the scattered and 
impressed electric fields involved in the equation, and neglecting cross polarized electric 
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currents (in view of the small cross section of the strip). The left hand side of (AO) 
represents the electric field at z=0 radiated by the electric current density j x . The Green’s 
function, g**, that appears in the kernel of this integral equation, is the electric field 
radiated by a short electric dipole lying at the interfile between the slab and the 
homogeneous half space for z<0, in the absence of the metallic strip. 

The electric current distribution is assumed to be characterized by a separable space- 
dependence with respect to x and y; i.e., 


j(x\y)=i(x) ■ j t (y) . (A1) 

The transverse y- dependence is chosen to verify the quasi-static edge-singularities 

Jt(y)~wTi (A 2) 

jt(y) possesses a closed- form Fourier transform, 

J t(ky)=FT{i t (y)}=J 0 (\wky), (A3) 

where J 0 is the Bessel function of zero order. The IE (1) is now transformed into the k x - 
Fourier domain as 


o°^ 

HK ) fGxxflCx, -y) j t (y) dy '= - E f (k x ), 

-oo 


(A4) 


where I( k x ) and G **(£,, -y) are the FT, with respect to the variable oii(x) and g^x, 
y), respectively. Equating the two spectra in A* it is apparent that 


Wx) D(t) ; (AS) 

where 

OO 

Dftx) =J G (k x , — y) j t (y)dy-j K f G(k x , k y )J Q ( \ w s k y )dk y (A 6) 

■°° -oo 

represent the spectral domain Green’s function of the microstrip line and G(k x ,k y ) is the FT 

of S(x, y) with respect to both the variables x and y. The electric current on the microstrip 
can then be expressed as: 




dk x 


(A 7) 
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